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ABSTRACT 



/A/'ms. The purpose of this paper is to assess the feasibility of axion detection by X-ray spectroscopy of the sun. 
Methods. We review the theory of axion-photon mode conversion with special attention to axions emitted in the 14.4 keV Ml decay 
of 57 Fe at the solar center. These then mode convert to photons in the outer layers of the solar envelope, and may in principle be 
detected subsequently as X-rays. 

Results. For axion masses above about 10~ 4 eV, resonant mode conversion at a layer where the axion mass matches the local electron 
plasma frequency is necessary. For axion masses above about 10 -2 eV, this mode conversion occurs too deep in the solar atmosphere 
for the resulting photon to escape the solar surface and be detected before Compton scattering obscures the line. At the (detectable) 
axion masses below this, the flux of mode converted photons predicted by axion models appears to be too low for detection to be 
feasible with current instrumentation. Nonresonant mode conversion for axion masses below 10~ 4 eV is also plausible, but with still 
lower predicted fluxes, since the axion coupling constant is related to it mass. 

Conclusions. Prospects for meaningful constraints on massive axion parameters from X-ray observations of this transition from 
the Sun do not appear to be promising. However parameters for massless counterparts (e.g. the "arion") may still result from such 
observations. It may mode convert in the outer layers of the solar atmosphere, but is not restricted by this to have a small coupling 
constant. 
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1. Introduction 

The basic form of the Lagrangian for the theory of strong in- 
teractions known as Quantum Chromodynamics (QCD) con- 
tains a pseudo-scalar term involving the gluon field strength 
tensor that violates symmetry under the actions of charge and 
parity reversal (CP violation). One practical consequence of 
this is the prediction of an electric dipole moment for the 
neutron in excess of the experimental limits by about ten or- 
ders of magnitude. The most generally accepted solution to 
this problem involves the introduction of a new global chiral 
symmetry (|Peccei & Quinn(1977a)| |Pecce i & Quin n(1977b)| l. 
Weinberg(1978) and Wilczek( 1978) pointed out that such a sym- 
metry, being spontaneously broken, necessitates the existence of 
a new boson, the axion, with mass in the range 10~ 12 eV to 
1 MeV, depending on the scale at which the global symmetry 
is broken. The axion mass is related to the symmetry breaking 
scale, f P Q, by (e.g. [Moriy ama( 1 993)) 



m = 0.62(l0 7 GeV// Pe ) eV. 



(1) 



Limits on the axion mass are placed by various astrophys- 
ical constraints. There are two most often discussed classes 
of axion; "hadronic" axions that do not couple to leptons 
and quarks at tree level (Ki m(1979)| |Shifman et al.(1980)] > 
and DFSZ axions that do ([Dine et al.(1981)| |Zhitnitskii(1980)|). 



hadronic axions may also have a mass of a few eV, bounded 
from below by treatments of axion trapping in SN 1987 A 
( |Burrows et al.(1990)) and from abo ve by studies of red gi- 
ant evolution ( Haxton & Lee(1991)| ). Studies of white dwarf 
cooling and pulsation also indicate DFSZ axion masses of or- 
der 10~ 3 eV (Tsern et al.(2010)). Cosmological arguments (e.g. 
|Kolb& Tur ner( 1994)) give a lower limit on the axion mass in 
the range 1Q~° - 10" 3 eV, with a mass of ~ 10 -5 eV supplying 
the measured cold dark matter density (Erken et a l.(201 1)) . 

The sun is the only astrophysical object besides SN 1987A 
from which neutrinos have been directly detected, and so it is 
natural to consider it also as a potential source of axions. While 
thermal processes at the solar center produce axions by Compton 
and Primakoff photoproduction processes, illustrated in Figure 
1, ( |Chanda et al.(198 8)), considerable attention has also been de- 
voted to a transition at 14.4 keV in the 57 Fe nucleus, which lies 
in a spectral region with essentially zero background from the 
nonflaring solar atmosphere. Being an Ml transition, it has a 
branching ratio for axion emission (DFSZ or hadronic) instead 
of photon radiation (Hax ton & Lee(1991)| l, and has a sufficiently 
low threshold to allow thermal excitation. An estimate of the ax- 



ion flux in this line has been given by Andriamonje et al.(2009) 
as Fa = 4.56 x I0 23 (g e /J) cm _2 s _1 , where g_J is the effec- 



Both types of axion are constrained to have masses less than 
about 10~ 3 - 10~ 2 eV based on analyses of the neutrino burst 
from SN 1987A ( |Burrows et al.(1989)| |Keil et al.(1997)) , while 



{sIn )~ cm-V 1 , where 
tive axion-nucleon coupling. Moriyama(1995) gives Fa — 5.3 x 
10 8 (m/leV) 2 cm~ 2 s _1 assuming the naive quark model (NQM) 
3.53 x 10~ 8 (m/leV). |Andriamonje et al.(2009)] 
10~ 8 (m/leV) (in the calculation of the 



where g e a f J 



give g'H = 1.48 x 
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Fe Ml axion branching ratio), which would give a flux of _ 1 I , ^ _ nf_\ + 1 ll e _i + «£.\ 2 + 4 
■ = 10 8 (m/leV) 2 cm" 2 s _1 but also discuss higher possible : 2 '^ 



SayyB 



(5) 



values for g e U . 






Fig. 1. Feynman diagrams illustrating the two Compton scatter- 
ing diagrams (left) and the Primakoff process (right). 



[Andriamonj e et al.(2009)| describe the CERN Axion Solar 
Telescope (CAST), designed to detect solar 57 Fe axions by ob- 
serving the 14.4 keV photons produced when the axions mode 
convert inside the detector. The parameters of CAST; a 9T mag- 
netic field extending over a length of s = 9.26 m, give a prod- 
uct Bs si 8 x 10 7 G cm, which is smaller than the product of a 
typical solar magnetic field and density scale height, and so the 
question of whether 14.4 keV photons might be detectable in the 
solar spectrum as a signature of in situ axion mode conversion 
becomes interesting to consider. 

2. Axion-Photon Interactions 

The axion-photon interaction Lagrangian is 



L = -^aF^F MV =g ay7 aE-B 



(2) 



where g ayy is the axion-photon coupling constant, a is the ax- 
ion field, F 1 ™ is the electromagnetic field tensor and F pv is its 
dual. The product F^F^ = -4E • B, which leads to the final 
form. Here, E is interpreted as the photon electric field and B the 
ambient magnetic field. Thus axions only interact with photons 
when propagating across the ambient magnetic field, and then 
only with the mode polarized along the magnetic field, known 
as the ordinary mode. 

The effective Lagrangian for an axion-photon system is then 

L=-E-e-E-^B-B + - (d p af - l -m 2 a 2 + gayy aE ■ B (3) 

where m is the axion mass and e = 1 - <y 2 /w 2 is the dielec- 
tric tensor in terms of the electron plasma frequency 10 p . We are 
working in the limit of small magnetic fields. By way of con- 
trast, Yoshimura(1988) includes the magnetic field terms in e. 
The equations of motion are 



e ■ E - §rE + g ayr aB = 
£2zM + (i_»4_ii) fl = o, 

cj- \ OJ- co- ) ' 

which for nontrivial solutions require 



(4) 



These two solutions represent axion and photon propagation 
across the magnetic field. We rewrite the photon-axion mass ma- 
trix from equation 4 in symmetric form as 



M = 



2(jJ 

gayyB 



gayyB 



2oj 



(6) 



Pre and post multiplying by a rotation matrix, 

cos 8 



cos 8 - sin I 
sin 8 cos 8 



M 



sin 8 

m" — 0)„ 



sin# 
cos 8 



^ ' cos 28 - g a yyB sin 28 ^ - sin 28 + g ayy B cos 28 



- sin 28 + gayyB cos 28 — - * cos 28 + g a7Y B sin 28 j 



7) 



which is diagonal for a mixing angle 6 given by tan 28 = 
2g aj yB(x>l [oj 2 p - ra 2 ). 

The evolution of the axion-photon mixed state with time t is 
described by 

\a (f)> = e- iE «' \a' (0)> cos 8 + e^' \y' (0)> sin 8 



= e- iE "' \a (0)) cos 2 8 + e- iE y< \a (0)> sin" 8 
+ cos 8 sin 8 (e""^ - e~ iE A \y (0)) . 

The probability of axion to photon conversion is then 



p { a -> y) = |(y|a (f))| = cos 2 8 sin 2 8 \e 



-iE„t 



sin 28 sin i — ^- L - 



which upon substituting for the mixing angle becomes 



P(a^y) = 



R 2 B 2 

oayy 



[E a - Ey) I 



gl n B^ + ((4-m^) IAa? 



(8) 



(9) 



(10) 



We put \E a - Ey} t = cjs (n a - n 7 j with n 2 y given by the two 
solutions in equation 5, and s the distance traveled through the 

medium. Then E a - E y ^ - w 2 ,) /a> 2 + \g 1 ayy B 1 j2 and 
P(a^y) =- f" r f 2 - 2 x 

sin 2 f ^{m 1 - cjff IAcj 2 + g 2 ayy B 2 
^ s 2 glyyB 2 /4 

(ID 

when (s/2) ^g 2 ly yB 2 + (w 2 - m 2 f /4w 2 ^ 0. 

The forgoing estimate assumed that the electron density in 
constant. When crossing a resonance, the rapid variation in 
tan 28 with u> 2 first greater than m 2 becoming equal and then less 
than m 2 , as the axion moves from high density to low density, 
means that the axion-photon oscillation cannot keep up with the 
changing Hamiltonian. In this limit, the equations of motion (4) 
admit solutions in terms of parabolic cylinder functions, yield- 
ing an axion transmission probabilit y of exp (-7ig 2 yy B 2 kh / '<^ p \ 
( |Cairns & Lashmore-Davies(1983)| |Cally(2006)} . The corre- 
sponding mode conversion probability is then 

P (a -» y) = 1 - exp (-7rg 2 ayy B 2 kh/cj 2 ,) - 7ig 2 yy B 2 kh/ oj 2 , (12) 

and as will be seen below, is generally larger than the nonreso- 
nant conversion probability given in equation 1 1 . 
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3. In Real Numbers ... 

We have used natural units, where % = c = 1, above. To 
put variables in units more familiar to solar physicists, we use 
|Andriamonje et al.(2009)l 



g ayy = 836xlO- 4 /fpQ GeV- 1 , 



(13) 



concentrating on DFSZ axions for the time being, where /pq is 
the symmetry breaking scale in GeV, and hence from equation 1 

m = 0.62(l0 7 GeV// Pe ) = 0.74(^ yr /10" 10 GeV- 1 )eV. (14) 

With these definitions 



g ayy B = lO-'^B^/lO-'OGeV- 1 ) 
= 1.3 X 10- 16 B(m/leV)cm-' 

where B is given in Gauss, and the conversion probability is 



(15) 



(16) 



p { a y) = 4.5 x l(T-"fi z (G) (m/leV) 2 sf. 

At a photon energy of 14.4 keV, 

a> 2 p /(o = 4.8 x 10- 2l n e cm- 1 
tn 2 /co = 3.53 WleV^cirT 1 . 

(17) 

For to > 10~ 4 eV and s > 10 8 cm (a typical solar atmosphere 
scale size) the small angle approximation for the sine in equa- 
tion 11 is not valid. Both terms in equation 17 are >> g 2 ayy B 2 
and so unless a> p to, the conversion probability is signifi- 
cantly reduced from the value in equation 16. The electron den- 
sity at which resonant mode conversion can occur is (see also 
|Yoshimura(1988)] l 



n e 7.35 x 10 20 (TO/leV) 2 . 



The electron scattering optical depth between this layer and the 
solar surface is 



n e o~ T h 



3 x 10 _1 V 



(19) 



where cr r is the Thomson cross section and h 500 km is the 
atmosphere density scale height. Figure 2 shows emergent line 
profiles at the top of the solar atmosphere after traveling through 
electron scattering optical depths of 0.1, 0.3, 1.0, 3.0, and 5.0, 
calculated following the treatment of Chandrasekhar( I960)] and 
Miinch(1948) forplane parallel atmospheres. As can be seen, the 
detectable line disappears once r ~ 5, meaning a maximum den- 
sity at which detectable mode conversion will occur is ~ 10 17 
cm -3 , with a maximum detectable axion mass of 1.2 x 10~ 2 eV. 
At higher opacities, a broad Compton scattered feature may still 
be seen. Figure 3 shows this feature at opacities of 5, 15, 50, 
and 150. It is visible as a 1-2 keV wide feature out to an opacity 
of between 15 and 50. This would increase the maximum de- 
tectable axion mass to about 4 - 7 x 10~ 2 eV, but would probably 
require more photons detected due to the broader feature. 

Equation 12 gives the mode conversion probability as P — 
7.5 x 10" 25 B 2 , independent of axion mass. For B ~ 10 3 G, P ~ 
10~ 18 , and the photon flux at earth is 10~'° Ow/leV) 2 cm~ 2 s 



-l 



assuming the axion flux calculated by Andriamonje et al.(2009) 
For axions with mass up to 10" 2 eV (the high end of the 
mass range allowed by SN 1987 A observations, and close 
to the limit allowed by Compton scattering in the solar at- 
mosphere), this leads to a count rate of order 10 _n s _1 
in a instrument of effective area 10 3 cm 2 (e.g. NUSTAR; 




0.00 28.80 57.60 86.40 
energy shift (eV) 



115.20 



Fig. 2. Line profiles of the emergent 14.4 keV 57 Fe line, origi- 
nally emitted as an axion, and mode converted at various optical 
depths in the solar atmosphere. The electron scattering opacities 
through which the various profiles have propagated are noted. 
The dashed curve gives the original unscattered profile. At an 
opacity of 5, the central narrow line has essentially disap- 
peared. 



o 
o. 



(18) £ 



o 
o 
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+ 
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0.00 57.60 115.20 172.80 230.40 
energy shift (eV) (Xio 1 ) 

Fig. 3. Same as Figure 1, but for electron scattering opacities of 
5, 15, 50, and 150. The r = 5 curve is truncated for numerical 
reasons. 



|Harrison et al.(2010)| |Harriso n et al.(2 013)) ; too low for detec- 
tion. Significantly higher conversion probabilities would be ex- 
pected in pulsar strength magnetic fields, (Yoshimura(1988); 
|Perna et al.(2012)) . 

For s ~ 10 s cm, the argument of the sine in equation 1 1 
becomes small for axion masses at the low end of the allowed 
range, to < 10~ 4 eV. In this case nonresonant conversion can 
occur with probability s 2 g 2 B 2 /4 =* 10~ 19 for B =* 10 3 G and 
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m 10 eV. This yields a flux at earth of mode converted pho- 
tons of 5 x 10 4 [g e / N f f . Taking g e / N f = 1.48 x 1(T 8 (>7i/leV) as 
before gives a negligible mode converted photon flux at earth 
of ~ 1CT 19 cm _2 s _1 with m = 10~ 4 eV. The lower axion mass 
leads to a significantly lower flux, through its coupling to g ayy . 

In hadronic axion models, the parameters g ayy and g e J^ are 
less tightly constrained. Yoshimura(1988) suggests g ayy may in- 
crease by a factor of ten, which would lead to a commensurate 
increase in the count rates estimated above, though still insuffi- 
cient to change the conclusions above. 



log(DFSZ axion mass (eV)) 

-6.13 -4.13 -2.13 -0.13 1.87 3.87^, 



o 
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-> o 
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DFSZ (res) 
Hadronic (nonres) 



DFSZ (nonres) 



s 

(Primakoff) 



.5? © 



and other constraints discussed therein. The dashed contour be- 
low this indicates the where this limit moves with the nondetec- 
tion of a count rate of 1(T 5 s" 1 by NUSTAR (a flux of 1(T 8 pho- 
tons cm _2 s _1 , and about 10 times the background count rate in 
this rate estimated from Table 2 in Harrison et al. 2013) for res- 
onant axion-photon mode conversion. The shaded contours in- 
dicate the allowed mass region (< 10~ 2 eV)in the lower left for 
DFSZ or hadronic axions. The lower dashed contour indicates 
the limit achieved by the same NUSTAR count rate for nonres- 
onant axion-photon mode conversion in the outer layers of the 
solar atmosphere. The rectangular boxes indicate the mass lim- 
its (< 10~ 4 eV) for this case. As can be seen, in either resonant 
or nonresonant cases, the expected flux at the appropriate axion 
mass is well below plausible detection limits. 

An exception might arise in the case of pseudo-Goldstone 
bosons arising from spontaneously broken family symmetry, 
which are naturally massless (Wilczek(1982) , but may oth- 



erwise have similar properties to the axion ( Anselm(1 988)i 



^ — ' |Berezhiani & Khlopov(1990a)| 



03 



10 O 



and|Berezhiani & Khlopov(199 



Berezhiani & Khlopov( 1990b) 



) give further discussion of the 



O 



-16.00 -14.00 -12.00 -10.00 -8.00 -6.00 



log(g ayy (GeV- 1 )) 



Fig. 4. The g ayy -g e JJ parameter space . The dark region at the top 
and right indicates the region ruled out by CAST (Andriamonje 
et al. 2009). The horizontal portion labeled (luminosity) comes 
from a constraint on the total Fe 57 axion luminosity as a frac- 
tion of the solar luminosity, and the vertical section labeled 
(Primakoff) comes from CAST constraints on axions produced 
by the Primakoff process, both from Andriamonje et al. (2009). 
The dashed contour below this indicates the where this limit 

by 



origin of such a particle, dubbed the "arion", in spontaneously 
broken family symmetry. Being massless, it may be detected 
following mode conversion anywhere in the solar atmosphere 
where the electron scattering opacity r < 1, which roughly cor- 
responds to the nonresonant case in Figure 4. In this case g ayy 

and g e JI are disconnected from the mass, and the forbidden re- 
gions of Figure 4. for DFSZ and Hadronic axions with mass 
above 10~ 2 eV no longer apply (the CAST forbidden region is 
still applicable). In such a case, NUSTAR might provide lim- 
its parameters competitive with those derived by other methods. 
For DFSZ axions, studies of the Primakoff process in massive 
stars limit g ayy < 0.8 x 10" 10 GeV" 1 flFriedland et al.(2012)) , a 
limit that would be accessible to solar observations. A similar 
limit comes from CAST observations of Primakoff solar axions 



(Andriamonje et al.(2007) i. 



NUSTAR (a flux of 10 photons cm s ) for resonant axion- 
photon mode conversion. The shaded contours indicate the "for- 
bidden" mass regions (> 10~ 2 eV) towards the upper right for 
DFSZ or hadronic axions corresponding to the values of g ayy 
and g e JJ (assuming the naive quark model). The lower dashed 
contour indicates the limit achieved by the same NUSTAR count 
rate for nonresonant axion-photon mode conversion in the outer 
layers of the solar atmosphere. The rectangular boxes indicate 
the "allowed" mass limits (< 10~ 4 eV) for this case. These "for- 
bidden" or "allowed" mass regions do not apply in the case of 
massless arions; only the CAST limits apply. Clearly, substan- 
tial parameters space is available to solar observations detecting 
nonresonant mode conversion of arions. 



These potential constraints are illustrated in Fig. 4., show- 
ing the g ajy - g e JJ parameter space. The top and right hand axes 
show the corresponding axion masses in the cases of DFSZ ax- 
ions (equation 1), and those coming from the naive quark model 
from g e Jj respectively. The dark region at the top and right indi- 
cates the region ruled out by CAST (Andriamonje et al. 2009), 



As mentioned above, axions may also be produced by ther- 
mal processes at the solar center (Chanda et al.(1988) l. Hadronic 
axions may also be produced by the Primakoff process. Compton 
scattering and electron-nucleus bremsstrahlung can also produce 
DFSZ axions. The Primakoff process alone produces around 10 
- 100 times more axions ( |Andriamonje et al.(2009)| l than the 
57 Fe decay, while the other two process increase the axion flux 
by approximately 10 3 (Moriyam a(1995)| |Derbin et al.(2011)) . 
The count rate so produced according to axion models would 
be ~ 10~ 8 s _1 at a mass of 10~ 2 eV, still too low for obser- 
vational feasibility, but possibly adequate to detect the arion. 
However these axions, and consequently the photons produced 
by their mode conversion, dominantly have energies in the 1 - 
3 keV range (Derbin et al.(2011)) . At these energies, the dom- 
inant opacity is photoelectric absorption rather than Compton 
scattering, and the depth of solar atmosphere from which mode 
converted photons may be detected is smaller. In this spectral 
band, the sun also emits a thermal bremsstrahlung spectrum, 
as well as spectral lines from active regions and flares, in con- 
trast to the 14.4 keV region, where no solar emission is ex- 
pected outside of flare. The fundamental problem would be that 
if a definitive detection is not expected, to be of value, the ob- 
servation needs to be designed so that a null result produces 
a further restriction of the observationally allowed parameter 
space. Hudson et al. (20 12) discuss progress and limitations in 
this approach. Davoudiasl & Huber(2006) suggest observing ax- 
ion mode conversion in the magnetic field of the earth, in such 
a manner that the earth itself blocks out the thermal continuum 
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from the sun. While this estimate uses the observationally al- 
lowed axion parameter space rather than a specific model to ar- 
rive at a count rate, it does have the advantage that a null result 
would be more readily interpreted. 

4. Conclusions 

The detection of photons from axion mode conversion within 
the solar atmosphere is practically limited to axion masses at the 
higher end of the mass window, ~ 10~ 5 - 10 -2 eV, (bounded 
by cosmological arguments from below, and by stellar luminosi- 
ties and SN 1987A from above). Detection of solar axions is also 
necessarily limited to masses below 10~ 2 eV due to the Compton 
scattering expected between the mode conversion layer (deeper 
in the solar atmosphere for higher mass axions) and the solar sur- 
face. This restriction means that current instrumentation is un- 
able to provide meaningful constraints on theoretical axion mod- 
els, because axion coupling constants scale as the axion mass, 
m, and so the mode converted signal scales as m 4 , leading to an 
intrinsically weak signal for low mass axions. A possible excep- 
tion to this might a massless counterpart to the axion, the arion, 
where g ajy and g"™ are now decoupled from the mass, and con- 
sequently the "forbidden" regions in Figure 4. do not apply. 
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